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Abstract 

Invariant tori are prominent features of symplectic and volume preserving maps. From the point of view of chaotic 
transport the most relevant tori are those that are barriers, and thus have codimension one. For an n-dimensional 
volume-preserving map, such tori are prevalent when the map is nearly "integrable," in the sense of having one action 
and « - 1 angle variables. As the map is perturbed, numerical studies show that the originally connected image of the 
frequency map acquires gaps due to resonances and domains of nonconvergence due to chaos. We present examples 
of a three-dimensional, generalized standard map for which there is a critical perturbation size, Sc, above which there 
are no tori. Numerical investigations to find the "last invariant torus" reveal some similarities to the behavior found by 
Greene near a critical invariant circle for area preserving maps: the crossing time through the newly destroyed torus 
appears to have a power law singularity at Sc, and the local phase space near the critical torus contains many high-order 
resonances. 
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1. Introduction 

Volume preserving maps are appropriate models for many systems including fluid flows [1, 2, 3, 4, 5, 6, 7], 
granular mixers [8], magnetic field line flows [9, 10, 11], and even the motion of comets perturbed by a planet on 
an elliptical orbit [12]. Volume-preserving dynamics has some similarities to symplectic dynamics; however, though 
every symplectic map is volume preserving, the converse is only true in two-dimensions. 

This paper is concerned with the effects of perturbation and resonance on invariant tori; such tori are especially 
common in the integrable case. In the context of Hamiltonian systems and symplectic maps, integrability is syn- 
onymous with Liouville's definition: a t/-degree of freedom system is integrable when it has a set of d, almost- 
everywhere independent, involutory invariants. The involution property of the invariants implies that they also gener- 
ate an Abelian group of symmetries that preserve the invariants, and therefore that the (compact) integral manifolds are 
tori. An integrable 2ii-dimensional symplectic map can be written (at least locally) in terms of angle-action variables 
(6, J) e T'"' X M'' as 

0' ^9 + Q(y) , 

(1) 

where we will take T'' = W' jl/' [13]. The concept of integrability is perhaps less well-formulated for the volume- 
preserving case. However, it seems quite natural to use Bogoyavlenskij's concept of broad integrability [14, 15]. 
Roughly speaking, a system of n ODEs is broadly integrable if it has k independent invariants and d commuting 
symmetries that preserve the invariants, where k + d - n. We propose a similar definition for a map — the integrable 
case corresponds to the form (1) as well, but now {9, 7) e T'' x [16]. 
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It is natural to study (1) on the universal cover, letting (x,z) e M"' x M*, so that (6, J) = (x mod l,z). We study a 
simple class of perturbations of (1), which for the lift becomes 



x' = X + Q(z') , 

(2) 

z = z - eg(x) , 

where the "force," g, is periodic, i.e., g(x + m) - g{x) for any m e 1/ . Since Q is evaluated at z' in (2), this map is a 
volume-preserving diffeomorphism for any smooth functions Q and gr It is exact volume preserving when g has zero 
average — or equivalently the zero Fourier component of g vanishes [17]. Since this is the only case for which (2) can 
have rotational tori — that is, tori homotopic to the tori of (1) — that are invariant, we will make this assumption. 
The rotation vector (or frequency) map: 

n : M* ^ K' (3) 

plays an especially important role in the dynamics of (1). For the integrable map, the forward orbit {(x,,z,) : f e N) of 
each initial condition (x(),zo) has a rotation vector 

w(x(),Z()) = Um ^^^—^ (4) 
r— >oo I 

given by f2(zo)- If ^(zo) is incommensurate, see §2, the orbit densely covers a cZ-dimensional torus; by contrast, when 
the rotation vector is resonant an orbit densely covers one or more lower dimensional tori. 

When the map (1) is perturbed, many of its li-tori are immediately destroyed; however, KAM theory implies that 
there will still be a large set of invariant tori if the perturbation is small enough and smooth enough and the frequency 
satisfies a nondegeneracy condition. This is rigorously true for k - d when the perturbed map is exact symplectic and 
satisfies a Holder condition (i.e. is C^^^ for some h^O), and Q satisfies a nondegeneracy condition such as the twist 
condition 

detDQ>c>0, (5) 

see, e.g., [18, 19]. 

Of course, when k < d, the number of actions is smaller than the number of angles and the matrix DO. is no longer 
square. The "nicest" case corresponds to rank(DQ) = k implying that the image of Q is an immersed A;-dimensional 
submanifold. 

Maps of the form (1) with k - 1, so-called one-action maps [20], have codimension-one invariant tori. KAM 
theory implies that codimension-one tori are robust features of nearly-integrable, analytic one-action maps [21, 22]. 
These theorems assume that O e C''^' and satisfies a nondegeneracy condition of the form 

det(Z)Q, D^n, . . . D'^Q) > c> , (6) 

similar to that used by Riissmann [23, 24]. 

By contrast, when I < k < d the invariant d - k dimensional tori of (1) need not be as robust. For example when 
there are two actions and one angle, almost all of the one-dimensional tori can apparently be immediately destroyed 
even under a smooth perturbation [25]. 

Though codimension-one tori are commonly observed in near-integrable, one-action maps, they are often destroyed 
by resonant bifurcations as the perturbation grows. The nature of these bifurcations is strongly influenced by the form 
of the frequency map. Even when the map satisfies (6) perturbations of (1) can have many of the features of symplectic 
maps that do not satisfy the twist condition [26], see §3. In §4, we will investigate which of these persistent tori is 
most robust. 



When d = k, (2) is symplectic with the form dx A dzj only if Dil and Dg are symmetric matrices. 
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2. Frequency Maps and Resonance 

A rotation vector w e M"' is "resonant" when there exists an (m, n) € Z'' x Z \ {0, 0) such that 

m ■ (jj - n . (7) 

The resonance module of a given w is a sub-lattice of 1/ defined by 

Ud) = {m e Z"' : m ■ w e Z) . (8) 

The dimension of this module (the number of independent ;«-vectors in X.) is the rank of the resonance. An incom- 
mensurate or nonresonant rotation vector corresponds to the rank-zero case, £, - [O]- The order of a resonance is the 
length of the smallest nonzero vector in X, though of course this depends upon the norm used; we will typically use 
the sup-norm. The set of resonant frequencies 



H = eW' : m ■ CO — n for some (m, n) e Z'' x Z \ {0, 0)} 



is the resonance web; it is a dense subset of W' . 

For the case at hand, d - 2 and k - \, the image of Q in (1) is a curve, as sketched in Fig. 1. As implied by 
KAM theory, we expect that, upon a general, smooth perturbation, only tori with "sufficiently" nonresonant rotation 
vectors (in the sense of a Diophantine condition) will persist when (1) is perturbed [21, 22], and that tori near low-order 
resonances will be destroyed most quickly. The set of Diophantine rotation vectors is a Cantor set, approximated by 
the white region in Fig. 1. The frequencies of the preserved set of tori correspond to a Cantor set near n(z): in contrast 
to the symplectic case, these need not be in the image of Q.. 

One way to find invariant tori is to numerically estimate the frequency map; orbits on which the frequency map is 
well-behaved should correspond to tori or periodic orbits [27]. Numerically, we can estimate the rotation vector (4) as 

(x>t{_xo,zo) = (9) 

for some T. The finite-time estimate (9) can be used to obtain an approximation of the image of the frequency map on 
the set of rotational tori. Since each rotational invariant torus will intersect every angle, we fix can the initial angles 
xq, and think of (9) primarily as a function of the initial action 

wr(JCo, ■) : ^ K''- 

The estimate (9) is rather crude since errors will decay as T"' even when the limit exists. There are more sophisticated 
methods for a single frequency [28, 29], but as these are based on continued fractions that rely on one-dimensional 
ordering, they cannot be used when d> \. Filtered Fourier methods [27] are also more accurate; however, these sufifer 
from the difficulty of identifying which peaks correspond to the "primary" rotation vectors. 



3. Tangency Normal Form 

GenericaUy the resonance web will intersect the curve Q(z) at a dense set of points. For a rank-one resonance the 
typical behavior is exemplified by a three-dimensional normal form obtained in [26] that is analogous to Chirikov's 
"standard map" [30] and Froeshle's four-dimensional symplectic map [31]. This normal form corresponds to (2) with 
{x, z) e X M' and the frequency map and force 

il(z)^{z + y-6+Pz^), ^^^^ 
g{x) - a sin(27rjici ) + b sin(27r.X2) + c sin(27r(.xi - X2)) . 

Note that £2 satisfies the nondegeneracy condition (6) since the vectors DD and D^Q. are never parallel: 

det(DQ, D^Q) = 2/3 . 
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Figure 1: Sketch of the resonance web, thickened by the Diophantine condition for a one-action map on T X R. The complement is the positive 
measure Cantor set of Diophantine rotation vectors. Also shown is the image, Cl{z), of a frequency map. 



Moreover, since g e C^(T^) and has zero average, KAM theory is applicable to (10). 

The main resonances for (10) occur when Q. is commensurate with one of the "forced" resonances in g\ the 
amplitudes a,h,c represent resonant forcing with m - (1,0), (0, 1), and ±(1,-1), respectively. To limit the number of 
parameters, we select 

y = i(V5- 1) ~ 0.61803 , je = 2, fl = Z7 = c=1.0, (11) 

and vary d and e for our computations. 

The locations of the forced resonances for the integrable case, t{m,n), are given in Table 1. Note that most 
resonances occur at two z-locations for each 6 because Q is parabolic. If Q. is tangent to a resonance line, there will be 
only one z-location, and the structure of the dynamics near this point resembles a nontwist area-preserving map [26]. 
Under the approximation that only one resonance is present and its amplitude is small, its width w — defined to be the 
maximum difference between the upper and lower separatrices — can be computed by standard techniques [26]. As an 
example, the phase portraits for 5 = 0.1 in Fig. 2 show four prominent resonances; the locations and widths of these 
forced resonances are given in Table 2. When e is small, as in the left pane, much of phase space away from these 
resonances is foliated by rotational tori; the region within the width of each forced resonance is predominantly foliated 
by "Ubrational" tori, the tubes in the figure. Of course, the dynamics near the resonance separatrices is chaotic, and as 
E increases, chaotic orbits become more prominent. 

A numerical computation of the approximate, perturbed frequency map using (9) for T - 10^ is shown in Fig. 3 
for the same parameters as Fig. 2 and with initial conditions along the line segment {(0,0, zo) '■ Zo e [-0.5, 0.5]). When 
E — 0.005, as in the left pane, the components of ojj follow those of Q(z) closely except near the forced resonances 
where the frequencies resonantly lock within a width w about z*. Locking intervals are visible for the (1,0, 1), (0, 1,0)^ 
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Figure 2: Some orbits of (2) witli (10) and (11), c5 = 0.1 ande = 0.005 (left), e = 0.015 (right). The outlined cube is [-0.5, 0.5]^ In the left panel, 
orbits trapped in the forced (1, 0, 1), (0, 1,0)^, and (1,-1,0)" resonances form gold, blue and grey tubes while the red and blue orbits are near the 
separatrix of the lower (0,1,0) resonance. In the right panel, the ( 1 , 0, 1 ) and (0,1,0)^ and the (0,1, 0)" and (1,-1,0)" resonance pairs overlap to 
form large chaotic layers. 
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Table 1: Locations of the forced resonances m ■ n(z*) = n for the frequency map (10), their full widths w, and values of 6 at which Sl(z) is tangent 
to the resonance. 

and (1, -1,0)" resonances; the fourth resonance (0, 1,0)" appears as a jump in loj instead of a locking interval because 
the chosen line of initial conditions passes near the hyperbolic invariant circle for that resonance, as can be seen in 
Fig. 2. Forced resonances in Table 2 with other values of n do not occur within the domain of the figure. 

Though the frequency components within a mode-locking interval appear constant on the scale of Fig. 3, they 
actually vary slowly, maintaining the resonance condition across the interval. 

Chaos is much more visible in the right panes of Figs. 2-3 where e - 0.015; in particular, the frequencies vary ir- 
regularly with zo for zo < -0.292, indicative of the strong chaos associated with overlap of the (0, 1 , 0)" and (1,-1, 0)" 
resonances. However, for zo e [-0.29,0.092] both components of loj vary relatively smoothly; the corresponding 
phase portrait shows many tori in this region. The values of each component in Fig. 3(b) are visibly shifted to the left 
from the unperturbed frequency (10); indeed the extremum occurs at zo = -0.03 so that (yj-(0, 0, -0.03) ~ f2(0). For 
Zo > 0.092, the frequencies alternately exhibit resonant trapping and chaotic layers. 

The chaos visible in Fig. 3(b) is less obvious in the image of the approximate frequency map z wj- shown in 
Fig. 4. The computed frequencies are close to the unperturbed parabolic curve (10); however, they fall visibly above 
Q(z) in the chaotic region near the overlapping (1,-1,0)" and (0,1,0)" resonances and they extend considerably 
beyond the left endpoint of the unperturbed curve, Q(-0.5). 

Computations of uj for four different initial angles are shown in Fig. 4. In most cases the four computed fre- 
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Figure 3: Numerically computed frequencies, ljj, with T = lO''. Parameters are the same as Fig. 2. The initial conditions are (0, 0, zo(7)) with with 
zqU) = ~0-5 + 500' ~ 1' ■ ■ ■ '500. The dots are the first (blue) and second (red) components of aij as a function of zo. ™d the curves are the 
components of f2(z) of (10). 
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Table 2: Comparison of the theoretical and experimental locations (estimated as the centers of the flat regions) and widths of the forced resonances 
for (5 = 0.1 and £ = 0.005 as in the left pane of Figs. 2 and 3. The width of the (0, 1,0)" resonance cannot be measured along this curve of initial 
conditions 

quency curves are nearly identical even though a trajectory with a given zo may have a significantly diff'erent computed 
frequency for each initial angle. Indeed, in the regular regions the computed frequencies for the same zq but different 
angles often differ by 10% even though they typically lie within 1% of the parabola (10). This is precisely as expected 
for a region foliated by rotational tori intersecting different angles at different heights. The difference between the 
initial angles is especially prominent near resonances. For example, the initial phases (0.5,0.5) and (0.3,0.2) cross 
the tube of librating orbits for the (0, 1,0)- resonance and they have mode-locked segments with ut ~ (|,0); how- 
ever the phases (0,0) and (0.5,0.5) appear to lie near the hyperbolic circle of this resonance and their corresponding 
frequencies jump as zo moves across the resonance. 

Variation of the image of the frequency map with s is shown in Fig. 5 for 5 = 0.0 and initial conditions on the 
line (0, 0, Zo)- For this value the parabola 0(z) is tangent to the (0, 1,0) resonance. When s = 0.001, the numerical 
frequency map (the black dots) is indistinguishable from the parabola Q(z), except for some small deviations when 
Zo 6 [-0.372, -0.344] where the orbits appear to be trapped near the (1, -1,0)- resonance with lot ~ ( jf, |f), when 
Zo e [0,0.08] where Q(z) is near the (0, 1,0) resonance, and when zo e [0.358,0.404] where it is near the (1,0, 1) 
resonance. Other resonance lines are shown in the figure, but when the perturbation is this small, they have little 
eff^ect on the image, indicating that most of phase space is foliated by invariant tori. For e - 0.01 much of the image 
still follows the parabola, but for s > 0.02 the image becomes increasingly irregular, correlating with regions of 
chaotic behavior in the phase space. Indeed, when s - 0.02 there appear to be no rotational invariant tori except for 
z e [0.091, 0.153]. When s - 0.03, there appear to be no rotational invariant tori for zo e [-0.5, 0.5]. We will explore 
this transition in the next section. 

Finally, the variation of the frequency map with 6 is shown Fig. 6. In this figure, gaps around the forced resonances 
are prominent, and smaller gaps around some of the higher order resonant lines, e.g. (1, 2, 0) and (2, 1,0) can also be 
seen. 
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Figure 4: Numerically computed frequency map for (10) with (11), e = 0.015 and S = 0.1, as in the right pane of Fig. 3. The four sets of points 
correspond to initial conditions on lines (xo,zo) with (xq) fixed as indicated and zo ranging from -0.5 to 0.5 in steps of 0.02. The dashed curve 
(purple) represents the image of unperturbed frequency map (10) for the same range of z. Resonance lines up to order 3 are shown in grey; the 
forced resonances are the thicker lines. 




Figure 5: Approximate images of the frequency map (9) for T = 10^ iterates for (10) with (1 1), 5 = and a = b = c = s. The value of e ranges 
from 0.001 to 0.03 as indicated. 
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Figure 6: Images of the approximate frequency map for (10) with (1 1), e = 0.01, T = lO', and varying S from -0.5 to 0.5 in steps of 0.1. The initial 
angle is (0, 0), and zo varies from -0.5 to 0.5. Note that the curve for 6 = -0.5 is equivalent to that for 5 = 0.5 under ai2 — * 1x12 + 1 

4. The Last Invariant Torus 

In 1968 John Greene developed a method for determining the "last" rotational invariant circle for an area-preserving 
map. Greene demonstrated that a circle with a given rotation number w exists when a sequence of periodic or- 
bits whose rotation numbers limit on to have bounded "residue". He provided compelling evidence that the last 
circle for Chirikov's standard map has the golden mean rotation number, that the "critical" invariant circle is not 
C^, and the phase space in its neighborhood exhibits self-similarity [32, 33]. This self-similarity leads to a renor- 
malization group explanation for the existence of invariant circles [34, 35]. More generally, it is hypothesized that 
the most robust invariant circles for twist maps have "noble" rotation numbers [36, 37, 38, 39]. Though there 
have been many attempts to generalize Greene's ideas to higher dimensions — in particular for symplectic maps 
[40, 41, 42, 43, 44, 45, 46, 47, 48, 49] — there is still no comparable theory that describes the robustness and de- 
struction of multi-dimensional tori. 

The question of which torus is the "last" is perhaps of less importance for symplectic maps because when d - 
k > 2, the tori do not form barriers to transport. When the tori have codimension-one, however, the question is more 
relevant. 

As £ increases in (2), rotational tori become less common. It appears that there is an Sc(J3, 6, y), such that whenever 
s > Sc there are no rotational invariant tori. Which is the "last" torus and what are its properties? 

To make this question easier to investigate, it helps to add a vertical translation symmetry to the map by assuming 
that there is a y e Z'' such that 

Q(z+l) = Q(z) + ;. (12) 
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In this case, (2) is isotopic to a three-dimensional version of the Dehn twist. For each orbit {{x,,z,) : f e Z) of the Hft 
of such a map, the set {{x, + k jt, z, + k) : t e Z} is also an orbit for any k e Z; consequently, the phase space structure 
of the map on T"^ x M is periodic in z- 
For our numerical experiments we set 

Q(z) ^(z + y,-6 + A sin^inz)) , (13) 

which has a parabolic form like that of (10) when |z| <sc 1. This map satisfies (12) with j - (1,0). Since the most 
robust tori of (10) appear to occur near z - Q, the last torus for (13) should be close to that for (10) for (3 - Ajp-. We 
choose g from (10), as before and set 

i = 1, y = i( V5 - 1), fl = ^7 = c = 1.0 , (14) 

so that ^2 varies over a unit interval. Periodicity in x implies that y and 6 can be restricted to any unit interval; thus 
we will study 6 e [-0.5, 0.5]. 




Figure 7: Frequency Maps for (2) with (13), witli (14) with T = 10^. The curves correspond to S ranging from -0.5 (black) to 0.5 (cyan) in steps of 
0.1. For the left pane, e = 0.01, and for the right, £ = 0.026. 



Frequency maps for (2) with (13) and the force g in (10) at two values of s are shown in Fig. 7. Since the map 
for e = 0.01 (left pane) is highly regular it indicates there are many invariant tori; however the irregularity of the 
frequencies for e - 0.026 (right pane), indicates that most tori have been destroyed. 

4.1. Crossing Time 

The simplest rigorous criterion that forbids the existence of invariant tori is the existence of "climbing orbits"; i.e., 
there are no rotational tori whenever there is an unbounded orbit. Since it is numerically impossible to verify that an 
orbit is unbounded, an implementation of this idea requires a bound on the vertical extent of rotational tori. Define the 
vertical extent of a set 7" c T'' X M as 

A(T') - max (z) - min (z) . 

{x,z)eT (x,z)er 

Rotational tori can be shown to have a bounded vertical extent for the case of symplectic maps with uniformly non- 
degenerate twist; indeed, "Birkhoff's Second Theorem" implies that a Lagrangian, rotational invariant torus on which 
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the dynamics is minimal must be a Lipschitz graph over the angle variables [50, 51, 52]. The vertical extent of any 
such torus is obtained from the Lipschitz constant. This theorem justifies "converse-KAM" theory and can be used 
both analytically and numerically to rule out the existence of rotational tori [53, 54]. Unfortunately, there can be no 
straightforward generalization of Birkhoff's theorem to a volume-preserving map of the form (2) since its rotational in- 
variant tori are not always graphs [26]. Such "meandering" tori (like those of nontwist area-preserving maps [55, 56]) 
often occur near resonant tangencies of the frequency map. 

Even though we know of no rigorous bound on the vertical extent of tori of (2), we will assume that one does exist, 
i.e., there is a A^ax such that A(7~) < A„,a i for all rotational invariant tori T. In this case, using the vertical periodicity of 
(13), rotational tori can be ruled out if there is an orbit that climbs a distance 1 +A„,ax- Numerical observations indicate 
that A„,a,i is quite small. The largest observed vertical extent, A(7~) ~ 0.2, occurs near a reconnection bifurcation [26]. 

Let the crossing time, tcr, of an initial condition be the minimal time over which z changes by 1 -i- Ai^ix- 

f„- = min{f > : k, - zol > 1 + Kiax) ■ 

We compute the crossing times with ls.max = 0.3 for a set of k - 10 initial conditions of the form (0,0, zo) obtained 
by varying zo in a narrow range about —y where there is a hyperbolic invariant circle (for e <k 1) due to the (1, 0, 0) 
resonance. To make the computation finite, the number of iterations for each initial condition was limited to at most 
lO'o. 

We observed that the crossing time grows rapidly as e decreases and appears to go to infinity at some critical value. 
Be, with a power-law singularity of the form 

f„ ~ Tie- K, s„ a) = ^ . (15) 
(e - Sc)" 

An example is shown in Fig. 8. If this form is correct, then Sc is the parameter at which the last invariant torus is 
apparently destroyed. 

When 5 = 0.1 and s > 0.0273 - ey, there were orbits in the trial set with f„- < 10'*\ thus Sc < su- However, when 
s e [0.02733,0.0277), ta- > 10'" for some orbits in the trial set; consequently, our computations cannot estimate the 
average crossing time in this range. Since the median can still be computed if the majority of trial orbits cross, we 
used the median crossing time to fit the data using ^ = 10 trials at each of N - 238 values of e e [0.02747, 0.1]. The 
spread in the data is estimated using the median absolute deviation 

tmad - medianllf ; - median(f,)l) . 

j=\...k ■' 1=1. ..it 

The three parameters, K, Sc, and a of (15) were fit by minimizing the squared log-error: 

" ^ Z Usi) - log Tier, K, e„ a)f . (16) 

The resulting fit, shown in Fig. 8 for 5 = 0.1, gives Sc = 0.02741, and an exponent a - 2.624 with an rms error 
E - 0.22. However, this estimated value of Sc is clearly too large since we observed crossing orbits at this parameter 
value (see also §4.2). A two-parameter fit, setting Sc = 0.0273, a more reasonable value, gives a - 2.772 and E - 0.24. 
The average crossing time also appears to fit the model (15), giving = 0.0274, and a - 2.617 for the full fit, and 
a - 2.753 for the fixed-Sc fit with similar rms errors. 

Though the data provide clear evidence for (15), it is difficult to compute the exponent accurately. In particular, 
some of the largest deviations from the fit occur near the singularity: it is possible that if computations could be 
done for tcr > 10"^, the exponent would be different. Nevertheless, it is intriguing that the observed exponent for the 
three-dimensional map is close to that for a critical noble invariant circle of a twist map (a = 3.05 [57]). 

Fits to the crossing time data for other values of 6 are similar and in each case there seems to be a power law 
singularity in the crossing time; however, the exponent varies considerably with 6. 
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Figure 8: Log-Log plot of the median crossing time (points) for 5 = 0.1 and A„,„, = 0.3 as a function of e - where = 0.0274L The error 
bars show the median absolute deviation. t,„ad, for 10 trials at each e value. The line is an unbiased fit to the logarithm of (15) that minimizes the 
log-mean square error (16). 

4.2. Graphical Exploration 

An alternative technique to estimate that also yields bounds on the location of rotational tori is to determine 
if orbits lie in a large, connected chaotic zone, or are trapped in a set with small vertical extent and thus presumably 
bounded by rotational tori. This method, primarily visual, uses the openGL-based iteration code to display the evolu- 
tion of an orbit, recall Fig. 2. This interactive program displays a trajectory as a continuously evolving "cloud" of its 
most recent 2'-' iterates. 

Using this method, it is quite easy to discover that the region containing rotational tori shrinks to a narrow range 
of initial conditions as e grows. For initial conditions on the line (0,0, zo) and e - 0.02, this range is zq e {zl,Zu) — 
(-0.0663, -0.0104) + 0.0001, as shown in Fig. 9. The experiments indicate that there are tori near the upper and lower 
bounds of this range. Orbits with zo e {zl,Zu) may be chaotic or lie on resonant, nonrotational tori; nevertheless, each 
such orbit appears to be trapped in {zl, Zu)- Moreover, the orbit of (0, 0, zl) is chaotic and was observed to move below 
z - -0.5, while the orbit of {Q,Q,zu) is chaotic and moves above z - 0.5. 

As e grows, the interval {zl,Zu) shrinks as shown in Fig. 9. From this data, we assert that there are no invariant tori 
for E - 0.0273 because we found a pair of orbits that "cross" any domain where tori could exist. In particular, the orbit 
of (0, 0, -0.056002) is chaotic and has min(z,) < -0.5, while the orbit beginning at (0, 0, -0.05612) has max(z,) > 0.5. 

Note, however, that the destruction of tori does not appear to be monotonic with e: in the figure the domain of 
tori shrinks almost to a point, (-0.053956, -0.053906), at £ = 0.0270, but then grows to (-0.056080, -0.053160) at 
e = 0.0271 before finally disappearing at e = 0.0273. Our observations indicate that there are no rotational tori for 
larger values of e, though we cannot rule out tori that might reappear for very small intervals of e. 

These experiments indicate again that e,. < 0.0273 for 5 = 0.1. For e = 0.02729, the orbit of (0,0,-0.0561), 
though weakly chaotic, appears to be trapped between rotational tori. This orbit provides an estimate of the rotation 
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Figure 9: The interval of rotational tori as a function of e for 5 = 0. 1 on the line (0, 0, zo)- There are no rotational tori for zo t (zl, Zu) for initial 
conditions (0,0, zo)- The vertical extent through xq = of some low-order resonances, labeled by (m\,m2,n), are also shown. The two (1,7,0) 
resonances appear to reconnect near e = 0.0255; unfortunately, their librational tori are destroyed before this happens and so the reconnection is 
hard to observe. 



vector of the "last rotational torus", from (9), 

oj * (0.619,-0.0859) 

It is difficult to compute this value more accurately without a more systematic method to locate tori; we plan to do this 
in a future paper [58]. 

Similar experiments for other values of 6 show that the critical e value ranges from about 0.0187 to 0.0347, with 
the largest observed value for 5 - -0.4. 

4.3. Frequency Maps 

As a final method to look for the last torus, we explore in more detail the numerical frequency maps shown in 
Fig. 7. Again, we focus on the case 5 -Q.l. 

It is easier to detect the irregularity of the frequency map from its components; these are shown for e - 0.015 
in Fig. 10 as a function of the initial z-value. Chaotic regions correspond to a scatter in toj values for nearby initial 
conditions. These are interspersed with intervals where coj appears to be more or less smooth indicating rotational or 
librational tori. 

Resonance gives rise to librational tori, recall Fig. 2. These are most easily seen when one of the components of 
(jj is constant, such as the (0, 1 , 0) and (1 , 0, 1) resonances labeled in the figure. However, even when both components 
of m nonzero — such as the (1,-1,0) resonance seen in the figure — the frequencies vary only slightly as the initial 
conditions cross a resonance so, on the scale of the lower pane of Fig. 10, these resonances also appear to be flat 
intervals. 

There are two intervals in the lower pane of Fig. 10 where coj appears to vary smoothly. However, this smooth 
variation of co is interspersed with resonant intervals, as can be seen in the upper pane, which is an enlargement around 
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Figure 10: Numerically computed frequency components (9) as a function of zo for £ = 0.015 and S = OA with xq = (0, 0) and 7 = 10^ (lower pane) 
and T = 10** (upper pane) for the frequency map (13). The relatively smooth variation of the frequencies near zo = -0.45 and -0.05 hints at the 
existence of rotational tori. The upper pane shows an enlargement of the dashed boxes around zo = 0. In the upper pane, the frequency components 
a>i and tJi are plotted on different scales as shown on the left and right edges of the graph. 



Zo — 0. Resonances can be systematically identified by looking for cases in which |ot ■ a> - n\ < c, for some small c. 
The figure identifies all resonances to order 12 that occur for more than one consecutive zo value, using c - 10 Note 
that the vertical scales of coi (shown on the left) and W2 (on the right) are different for the upper pane. On the scale 
of this figure, the component wi is still appears to be constant across a resonance; however, the variation of W2 is now 
more visible since it is displayed on a smaller scale. Since (l)2 varies and m ■ w = « in the resonance, wi must vary as 
well. 

Though there were two small intervals in Fig. 10 containing invariant tori. Fig. 11 suggests that when s = 0.02 
there are rotational tori only in a small interval around z = -0.05. The upper pane is an enlargement of this region. 
Recall from Fig. 9, there are no tori outside the interval (zl,Zu) = (-0.0663,-0.0104); this range is also indicated 
in the figure by the dashed vertical lines. The endpoints of the region of smooth variation of oj clearly correlate with 
the interval {zl,Zu), confirming that smooth variation of oj indicates the existence of tori. As noted earlier, there are 
chaotic orbits trapped by tori within {zl,Zu)- 

For e - 0.025, (zl, Zu) - (-0.054260, -0.045980), as shown in Fig. 12. It is noteworthy that much of this interval 
is filled by resonant or chaotic orbits; this can be seen in the uppermost pane of this figure. The Cantor set of tori for 
this value of e is sparse indeed. 

When E - 0.026, there are tori only in the interval (zl, Zu) - (-0.055500, -0.047900). Unlike the previous figures, 
the upper bound that we found in §4.2 is not confirmed by the frequency map in Fig. 13: the frequencies near this 
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Figure 1 1 : Numerically computed frequencies, CJ7-, as a function of for £ = 0.02 and (5 = 0.1 with xq = (0, 0) and T = lO' (lower pane) and 
T = 10* (upper pane) for the frequency map (13). There appear to be invariant tori only in the range zo £ (-0.0663, -0.0104), and the upper plot is 
an enlargement of this region. Resonances up to order 12 are labeled. 



point vary irregularly, at least on the scale of the figure. It is possible; however, that there is a nearly isolated invariant 
torus at this location. 

When e is increased to 0.027 or beyond, it is diflicult to compute a meaningful enlargement of the frequency map 
for {zl,Zu) as this requires increasing T beyond 10^. Nevertheless, the frequency maps confirm the conclusions of the 
previous methods that ~ 0.0273. 



5. Discussion 

We have shown that the break-up of two-dimensional, rotational tori in a three-dimensional, volume preserving 
map has a number of features similar to those found for the area-preserving case. Namely, the "last torus" appears 
to be the limit of a family of resonances. However, we have not been able to discover if there is an analogue of the 
dynamical self-similarity discovered by John Greene for the standard map. We also saw that the time for an orbit to 
cross a region has a power-law singularity just hke that discovered by Boris Chirikov for the two-dimensional case. 
The implication is that destroyed invariant tori may leave remnants like the cantori discovered by Ian Percival for 
symplectic maps. However, as of yet, we know of no general construction or even any examples of such invariant 
cantori sets for the volume-preserving case. 

There are many open questions. Birkhofi"s theorem implies that the set twist maps with invariant circles is closed 
since the circles are Lipschitz [59, §5.4]. Is this true for the volume-preserving case? Can one bound the vertical 
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Figure 12: Numerically computed frequencies, ojt, as a function of zo for c = 0.025 and 5 = 0.1 with xq = (0, 0) and T = 10^ (lower pane) T = lO' 
(middle pane) and T = lO' (upper pane) for the frequency map (13). There are no invariant tori outside the interval z = [-0.05426, -0.04598] 



extent, A(7~), of rotational tori? Is the set of parameter values for where there are rotational tori itself closed? Is 
there an analogue to the robustness of circles with noble rotation numbers to the higher-dimensional case? Since the 
noble numbers are quadratic irrationals, it has long been speculated that cubic irrationals could fill this role in the two- 
frequency case. As Mather demonstrated, it is easy to show that the standard map has no rotational invariant circles 
for A; > I [53]; can one find an explicit ball in parameter space outside of which the family (2) has no rotational tori? 

There are two natural methods for numerical approximation of invariant tori. The first is to look for sequences of 
periodic orbits whose rotation vectors limit to one for a given torus. The second is to use a Fourier series approximation 
for the embedding : T'' — > T'' x M*^ that conjugates the dynamics on the torus to the rigid rotation t-^ 6 + a>. These 
complementary methods should lead to more insight into the break-up of tori, and we hope to report on them in the 
future [58]. 
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